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■ We compute two infinite series of tree- level amplitudes with a massive scalar pair and an arbitrary 
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number of gluons. We provide results for amplitudes where all gluons have identical helicity, and 
amplitudes with one gluon of opposite helicity. These amplitudes are useful for unitarity-based 
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K*" ' one-loop calculations in nonsupersymmetric gauge theories generally, and QCD in particular. 
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I. INTRODUCTION 



Explicit computations of Standard Model processes will play an essential role in probing 
beyond it. The short-distance experimental environment at hadron colliders requires the 
computation of many processes involving QCD interactions. Tree-level computations do not 
suffice for these purposes. The coupling in QCD is sufficiently large, and varies sufficiently 
with scale, that even a basic quantitative understanding [1] requires the computation of 
next-to-leading order corrections [2] to cross sections. One-loop amplitudes are of course an 
essential ingredient of such corrections. 

Within the unitarity-based method [3, 4, 5], one can decompose one-loop QCD amplitudes 
into contributions corresponding to A/" = 4, A/" = 1, and remaining contributions. For 
amplitudes where all external particles are gluons, this decomposition for the gluon-loop 
contributions of color-ordered amplitudes takes the following form, 

A„ = <=^-4<=i+Ar'^^ (1) 

That is, the remaining contributions correspond to scalars circulating in the loop. (The 
contributions of quarks circulating in the loop can also be written in terms of A'^^^ and 
^scalar The supcrsymmctric contributions can be computed by performing the cut algebra 
strictly in four dimensions, with only the loop integrations computed in D = 4 — 2e dimen- 
sions. The 'scalar' contributions require that the cut algebra, and the corresponding tree 
amplitudes fed into the unitarity machinery, also be computed in D dimensions [4, 6, 7, 8]. 
At one loop, computing a scalar loop in D dimensions is equivalent to computing a massive 
scalar loop, and then integrating over the mass with an appropriate weighting. 

The computation of tree-level amplitudes with massive scalars is thus of use in the unitar- 
ity method for computing massless loop amplitudes in nonsupersymmetric gauge theories. 
The simplest such amplitudes, with up to four gluons of positive helicity, were computed 
by Bern, Dixon, Dunbar, and one of the authors [8]. Recently, Badger, Glover, Khoze 
and Svrcek (BGKS) have given [9] a set of on-shell recursion relations for amplitudes with 
massive scalars. They used the recursion relations to compute all massive scalar amplitudes 
with up to four external gluons. These relations extend the tree-level on-shell recurrence 
relations of Britto, Cachazo, Feng, and Witten [10, 11]. On-shell recursion relations have 
also been applied to tree-level amplitudes by Luo and Wen [12] and by Britto, Feng, Roiban, 
Spradlin and Volovich [13]; to tree- level gravitational amplitudes by Bedford, Brandhuber, 
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Spence, and Travaglini [14], and by Cachazo and Svrcek [15]; and to massive vector and 
quark amplitudes by Badger, Glover, and Khoze [16]. The concept has also been applied 
to certain loop amplitudes by Bern, Dixon, and one of the authors [17]; and to the direct 
calculation of some integral coefficients by Bern, Bjerrum-Bohr, Dunbar, and Ita [18]. These 
relations grew out of investigations [19, 20, 21, 22, 23, 24] motivated by Witten's topological 
twistor-string description [25] (as a weak-weak dual) of the Af — A supersymmetric gauge 
theory. The roots of this duality lie in Nair's description [26] of the simplest gauge-theory 
scattering amplitudes in terms of projective-space correlators. For use as building blocks in 
loop amplitudes, we need analytic expressions for the massive-scalar amplitudes. While the 
original tree-level recursive approach [27, 28, 29] pioneered by Berends and Giele is more 
efficient for purely numerical purposes, the on-shell approach seems better suited to pro- 
viding analytic formulae. (We discuss the computational complexity of the Berends-Giele 
recursion relations in the appendix.) Bern, Dixon, and one of the authors have recently 
presented [30] a 'unitarity-bootstrap' approach which combines four-dimensional unitarity 
cuts with use of a recursion relation for the rational terms. The amplitudes we compute 
here should be useful for proving the factorization properties in complex momenta required 
for the recursion relations part of this approach. 

In this paper, we will provide formulae for two classes of amplitudes with two color- 
adjacent massive scalars and n gluons, where the gluons all share the same helicity, 
or where one gluon has opposite helicity. In the next section, we document the nota- 
tion and conventions we use; in section III, we derive a form for the all-plus amplitude, 
An{ls, 2"*", . . . ,{n — 1)"*", Us). In section IV, we derive an expression for the amplitude with 
a negative-helicity gluon adjacent to a scalar leg, An{ls, 2+, . . . , (n — 1)~, Ug), and extend it 
to one negative-helicity gluon in an arbitrary position in section V, followed by concluding 
remarks. 

II. NOTATION AND ON-SHELL RECURSION RELATION 

We will write our expressions for various amplitudes using the spinor-helicity formalism. 
The formalism makes use of spinor products. We follow the conventions of the standard 
QCD literature, so that 

(jl) = 0-|/+) = u-{kj)u4ki) , Ul]^ {j^\n = u4k,)u4ki) , (2) 
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where u±{k) is a massless Weyl spinor with momentum k and positive or negative chiral- 
ity [31, 32]. Wc take all legs to be outgoing. The two spinor products are related, with 
[ij] — sign(A;°A;°) {j i)* so that, 

{ij)[ji\^2ki-kj. (3) 

(Note that the bracket product [ij] differs by an overall sign from that commonly used in 
twistor-space studies [25] and also in ref. [9].) 

In the amplitudes we consider, we will also encounter sums of momenta, 

Ki...j = ki + --- + kj, (4) 

but in this paper the indices should not be interpreted in a cyclic manner; if i > j, Ki...j = 0. 
Let us also define a notation for a sequential product of spinor products, 

{{Ji ■ ■■32)) = (Ji (Ji + 1)) ((Ji + 1) (Ji + 2)) • ■ • - 1)3,) , (5) 

and one for the mass-subtracted square of momentum sums, 

j j 

Li-j = ka-h, (6) 

a=i b=i,^a 

SO that if k\ = m^, and k,, ■ ■ ■ , kj are massless, for example, then Li...j = Kf...j — m^. 

The on-shell recursion relations make use of complex momenta, obtained by shifting 
spinors corresponding to massless momenta. (One can also shift momenta of massive parti- 
cles, but we will not need to do so in this paper.) A {j, I) shift is defined by, 

Ij-) ^ Ij-) - z\n , 

lo ^r)+z\r), (7) 

with the remaining spinors unchanged. It gives the following shift of the momenta, 

k^ -^k,^k^i^)-k^-'^{r\r\i-), 

k^ ^knz)^k'; + '-{r\7'\i-). (8) 

To obtain an explicit expression using the on-shell recursion relations, we must choose a 
definite shift. The relations then express an amplitude in terms of a sum over all inequivalent 
contiguous partitions of the external momenta into two sets, each of which contain exactly 
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one of the shift momenta. These partitions can thus be thought of as corresponding to sums 
of cychcly-consecutive momenta. There is also a sum over the hehcities of the intermediate 
leg. The basic relation [9, 10, 11] is, 

partitions P h=± 

X ' ^ Anikp^, ...,ki,...,kp_^, P-""). (9) 

± TTlp 

In this equation. Pi stands for the first momentum in P, P_i for the last momentum in P, 
and P for the complement of P (all the remaining momenta). The scalar legs (always the 
first and last in this paper) will be indicated by an 's' subscript. The mass of the particle 
in the factorized channel is mp, and P is given by momentum conservation, 

p = p + 5k^{z) = kp, + --- + kp_, - ^-{r\r |r>. (10) 

In each channel, a different value of z is used here, and in eqs. (7) and (8). It is given by 
the on-shell condition P^ — mp = 0, 



(nrr) 

As discussed in refs. [9, 11], we must ensure that A{z) vanishes at large z for the recursion 
relations to be valid. The choices we will use have the required property, as was shown in 
these references. 

III. THE A{ls, 2+, . . . , (n - 1)+, n^) AMPLITUDE 

In this section, we provide a result for the amplitude with an arbitrary number of positive- 
helicity gluons. The amplitude with all negative-helicity gluons can be obtained by spinor 
conjugation. Let us begin at the end, by writing down the answer, 

i 



^n(ls,2 , . . . , (n 1) ,ns) — ,, , 

Ln/2J-1 n-3 j-1 I -\ tJ- i _\ 



'<'i=Wi_l+2 



In this equation, = k\ = k\ is the mass squared of the scalar, and \_x\ denotes the largest 
integer smaller than or equal to x. 



The multiple sums in eq. (12) can be written less succinctly and perhaps less forbiddingly 



as, 



n-3 

E 



wo=L, 
W4 =n—l 



n—3 n—3 

= E E ■ 

wi=3 W2=wi+2 



n—3 

E 

■U),_l='U),_2+2 



(13) 



For example, in the four-point case, the product over r is absent, and we obtain, 

im^ [3 2] (2- 11^1 1 2") im^ [32] 



Li2 (2 3) L34 



L12 (2 3) ' 



in the five-point case. 



L12 (2 3) (34) L45 

and in the six-point case, where there is only one non-trivial w variable, 

' ; (S"*" I 1^41^3 1 2 )" 



zm: 



5+|^2..Ji|2-> 



(14) 



(15) 



(16) 



Li2 (2 3)(34)(4 5)L56V^~' i.-.-.--m- / ^ 

all in agreement with refs. [8, 9] (up to an overall phase for the second reference). 

Having written down the above ansatz, we now prove it using the on-shell recursion 
relations as given by BGKS. For this purpose, let us choose a (3, 2) shift using eqs. (7) 
and (8). (That is, choose the 'reference' momenta to be ks and k2-) 

In eq. (12), the limits on the sums ensure that i < j in sums of momenta Ki...j. We 
can remove those limits, so long as we take Ki...j = ii i > j (rather than interpreting the 
indices in a cyclic sense), and we shall do so. We can then rewrite eq. (12) in a form which 
will be more useful for the proof. 



^„(l„2+,...,(n-l)+,n, 

L«/2J-1 n-3 j-l 

X E E (n 



Li2((2---(n-l)))L(„_i)„ 



{^i} 



Ll...(wr-l)Ll---V 



(17) 

Wi=n—l 



We will proceed inductively, assuming that the ansatz (12) holds for the {n — l)-point 
amphtude (n > 4), and showing that it holds for the n-point one. The on-shell recur- 
sion relations [9, 10, 11] tell us that the n-point amplitude can be written as a sum over 
all factorizations, with the two shifted momenta attached to different amplitudes in each 
factorization. 

No factorization can isolate k2 in a purely-gluonic amplitude. In addition, any factoriza- 
tion that isolates k^ in a purely-gluonic amplitude leads to a vanishing contribution because 



the amplitude vanishes. For factorizations where is isolated in a four- or higher-point 
amplitude, this is immediate, because all such amplitudes vanish, 

^,+i(i^3l^,3+,...,j+) = 0. (18) 

In the three-point case, the factorization is in the [3 4] channel, but here the relevant ampli- 
tude — As{Kf4^, 3"*", 4+) — also vanishes. 

This leaves us only with factorization in which both amplitudes have massive scalar 
legs, and therefore in which the factorized leg is a massive scalar. There is only one such 
contribution, so what we are seeking to prove is that, 

Anih, 2+,...,{n- 1)+, Us) = ^3(1., 2+, -Ki2s)-^2-^^An-i{ki2s, 3+, . . . , (n - 1)+, Us). 

(19) 

The three-point amplitude was given by BGKS [9] , 

A3(l.,2+,3.)=zi^:l|p. (20) 

Its form depends on an arbitrary reference momentum q, but its value is nonetheless inde- 
pendent of it on shell. 

Using this expression with q — k^, and rewriting denominators using momentum conser- 
vation, we have for our starting point. 



L12L123 ((2 ■ ■ ■ (n- 1))) L(„_i)„ 

L(n-1)/2J-1 „-3 j-1 

X E (-™;)' E (n 



{Wr I k{wr+l)-iwr+i-l) \'^r+l 



^ {«'i}i=i =4 



X (^/;i+|^3...K-i)^i2^3l^i|2-> 



, (21) 



where we have also used the fact that |2 ) = |2 ) and (3 | = (3 |. Now, 



±x2hW^-) = 2X12 •A;3 1^1 |2-)-^3i'i2l^i|2-) 

= U2zh\2-)-m%\r), (22) 



thanks to the on-shell conditions in the three-point amplitude. 
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Writing out the two terms, we have 



Li2((2.-.(n-l)))L(„_i) 

L{n-1)/2J-1 n-3 J-l 



E (-<r E (n 

L(n-1)/2J n-3 j-l 

+ E E (n 



Wr I {■Wr + l)-{Wr + l-l) l'"^r+l' 



-) 



3---(«;i-l)i 



^i|2-> 



(23) 

Wi=n—1 



We can rewrite the last factor in brackets, 

(W2+| t3-{w2-l)h |2") (^'^rl ^4-..(^2-l) k2") 



L4... 



'4-n 



tu,=n— 1 



> . 



L 



4---n 



<^i+|^2..>i-i)^i|2"> 



(24) 



«)i=3 



We also note that because of the on-shell conditions on the three-point amphtude, 
2^2 ■ ki — — (2~| |2~), so that we can replace (u^i+l ^^-(wi-i)^! |2~) in the last factor 

by («;i+|^2...K-i)l^a|2-). 

Furthermore, we can extend all but the Wi sums in eq. (23) down to w = 3, since the 
summands will necessarily vanish. The second term will then supply the Wi = 3 terms, with 
some left-over pieces. 



Li2((2---(n-l)))L 



(n—l)n 



L(n-1)/2J-1 



n-3 j-l 



Wr+l)---{Wr+l-l) 



E (-^'.y E n L L, . 

■> {«'iK=i=3 



2---(«)i-l)i 



^i|2-> 



io,=n— 1 



H-<v E (n 



n— 3 



1) Fr+1 



-)(«^i+|^2...K-i)l^i|2-> 



j=L(n-l)/2J 



(25) 



The last term is absent for n > 4 (because the constraints are incompatible). For the 
penultimate term, we consider the even and odd cases separately. 
If n is even, 

L(n-1)/2J-1 \n/2\-2 

E = E . (26) 
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and the penultimate term completes the sum in the first term to up to j = [n/2\ — 1. Note 
that in the penultimate term, Wi can in any event take only the value 3, because while we 
have Wr > 2r + 1 , the presence of n/2 — 2 different summation indices Wr with Wr > Wr-i + 2 
requires that Wr take exactly the value 2r + 1. 

If n is odd, on the other hand, then the penultimate term in eq. (25) vanishes, because 
we have (n — 3)/2 different summation indices Wr, and the last one would have to obey the 
incompatible constraints w^(n-3)/2 < n — 3 and W(^n-3)/2 > n — 2. 

We thus obtain eq. (17), as desired. We have verified this equation numerically, using a 
set of light-cone recursion relations of the conventional kind [28], through n — 12. 

IV. THE A{ls,2+,...,{n-2)+,{n-iy,ns) AMPLITUDE 

We turn next to an amplitude with one negative helicity, adjacent to a scalar leg. There 
are two such color-ordered amplitudes, related by reflection. 
Here, we shall use an (n — 1, n — 2) shift, 

l(n-l)-) -|(n-l)-)-z|(n-2)-), 

|(n-2)+) ^|(n-2)+) + ^|(n-l)+), (27) 
The recursion for the target amplitude now has two terms, 

An-l{ls, 2+,..., (^)\ ' _m^ M-Kl_,^ns^ (^)-, n,) 

+A_i(l„ 2+,..., i^(t_3)(n-2)' (^)-,ns) (28) 

X M-K^n-S)in-2r ' 3) + , (^)+)- 

-^(n-3)(n-2) 

We can evaluate the first term directly, 

Ti(i„ 2+, . . . , (^r^)+, (;r^)-, n,) = 

An-l{ls, 2+, . . . , (^) + , Vi).J ^. ' _^2 M-Kl_,^ns^ (^)-, n,) 

-^(n-l)n '"'s 

^ i ((n- 1)- 1 \q-) ^^g^ 

L,, ((1 ... (n - 3))) {{n - 3) (r^)) L^^^^^_^^L~^ [q {^l)] 
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L(n-1)/2J-1 „-4 j-1 

X E E (n 



Ll...(wr-l)Ll---Wr 



-) (^i+|^2...K-i)^i|2-> 



■Wi=n—1 



Li2 ((2 • • • (n - 3))) ((n - 3) | ^(„_2){n-i)4^ (n-2)(n-l) l(^ ~ 2) ) L(„_2)...n-^(n-l)r 

L(n-1)/2J-1 n-4 J-l , „ _, , „ _\ 



- E 



-m 



2V 



E n 



For the second term, 



Ll...(wr-l)Lh 



-) <^i+|^2..>,-i)^i|2-> 



T2(l„2+,...,(n-2)+,(n-l)-,n,) = 



(30) 



X -^2 M-Kln-3)in-2)^ " 3) + , - 2)" 

-^(n-3)(n-2) 



we can iterate the recursion relations; suppressing all arguments but the number of legs, we 
obtain the following structure. 



(31) 



= + ^ [Ti(n - 1) + ^ [T,(n - 2) + ^^^.3]] 



In this iteration, the K momentum at any stage will be one of the shifted momenta at 
the following stage; to make this explicit, it will be helpful to define. 



— 

[1] ~ '^ra-\i 
-^S] ^ -^m-j,^y_i] = Kn-3 + 



(/c^_j + Ky_i])^ (m |7^l^b-il ) 



2(m \i^_^\K^-x^ ) 



(32) 



(/c„_,_i + %])2(m-|7^|ir[,f) 



2 (m I iyn-j-\ |%] ) 

so that -^[j-i] + — -^[j-i] = 0. For the specific configuration considered in this section 
m = n — 1, and hence we have 



^[1] 



(/c.-,-i + 4-1])^ ((n - 1)-| 7^ 

2((n-l)-||^„_,._i|i^[,_if) 



(33) 



10 



b] b] 



{K_,.2 + Ki,])mn-l)-\r\Kun 



2((n-l)-|4-,-2l^br) 
Putting the arguments back into eq. (31), we obtain an explicit expression for An, 

^„(l„2+,...,(n-2)+,(n-l)-,n,) = 

^^ a3(-^t;p(n-i-r)+,f+_,]) 

j=l r=2 (n-2-r)-{n-2) 

xTi(l„ 2+, . . . , (n - 2 - jy, k+, (^T^)-, n,). (34) 

To proceed, we must simplify the product of three-point amplitudes. We can first rewrite 
a3(-^p],(ri-3)+,fj]) _ 



-^(n-3)(n-2) 



[(n-3)%]r 



[(-i^p]) (n - 3)] [(n - 3) K^,]] [K^-,] (-i^pj)] Kl_,^^^_ 



■2) 



((n-i)(-;^p])r 



((n - 1) (n - 3)) ((n - 3) (n - 2)) ((n - 2) (n - 1)) 
and then by induction, we can show that 



(35) 



r=2 



{m-r)-{m-l) 



{m (m — j + 1)) (((m — j + 1) ■ ■ ■ m)) 
where for the case considered in this section m — n — 1. 

Using this result, we obtain our final formula for the amplitude, 

yl(l„2+,...,(n-2) + ,(n-l)-,n,) = 

z((n-l)-|^i4|(n-l)+)' 



(36) 



((2 • • • (n - 1))) Kl^^^_,) {{n - 1)1 ^2...(n-i)^i |2+) 

i{{n-l-2) (n-l-l)) 

+ ^ ((2 • • • (n - 1))) Li2i^(l_,_i)...(„_i)i^(„_,_i)...„ 

{{n - 1)-| $:^n-i-i)-{n-i)l^n \{n - 1)" 



(37) 



X 



{{n-l-2) I ^(n-i-l)...(n-l)^n l(^ " 1)"^) ((^ " ^ " 1) I ^ {n-l-l)-in-l)L I (^ " 1)" 

L(n-0/2j-l n-l-3 J-1 , , , _^ 



-m 



i: (n 



Wi=Wi_i+2 



Ll...(wr-l)Ll---Wr 

X (^i^i+|^2..>i-i)^i|2"> 
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where we have separated out the m-independent terms. The reader can readily verify that 
these reproduce the required maximally helicity-violating (MHV) amplitude in the massless 
limit. 

We can rewrite eq. (37) in a slightly more compact form as 

A(l„2+,...,(n-2)+,(n-l)-,n,) = 

i{{n-l)-\Mn\{n-iyf 



112... {n- 1))) Kl,^^_,^ {{n - 1)1 ^2...(n-i)^i |2+) 

n-4 



i((n-l-2){n-l-l)) 
+ Z. ((2 . . . (n - 1))) i^i2i^(l_,_,)...(„_,) ^ ' 

^ {{n - 1)- 1 ^(n-;-i)...(n-i)l^n l(^ " 1)+)' G{2, n-l-l;n-l) 

""{{n-l- 2)-| ^(„_,_i)...(„_i)|^J(n - 1)+) {{n-l- 1)-| |(n - 1)+) 

by introducing 

L(6l-ai+l)/2j bi-2 i-1 / _, jrf^ I _\ 

G{ai,bi;m) = {-mj 2^ [[[ '-^ — '- 1 

7^ ,^-1 V=i / 



X (wi+l |ai" 

= *^ 1 — 1 ; = & 1 

In the four-point case, only the first term in eq. (37) is present, and we obtain, 

I 2 



(^+l = <m-|^f(,^...^, 



(39) 



^4(l.,2+,3-,4,) 



\23)Kl,{3-\Uim 

-i'^W^: (40) 

-"^23-^12 



in the five-point case. 



Kl^, (23) (34) (2-||^i^23|4+) -^12-^45 

[34] (4-||^5^34|2+)' 

and in the six-point case, 

^(5-| ^-1^234 15+)' ^ ^m2(5-|^6|4-)'(4+|jt23^i|2-) 



(41) 



(2 3) (3 4) (4 5) (2- 1 ^1^234 |5+) ^^12^456^56 (2 3) [4 5] (5"! ^^45 |3 
imj (2 3) (5-|fe^34|5+)'(5-|^34|2-) 
Kl,,L,, (2 3) (3 4) (4 5) (5"! ^^345 |2+) (51 |3+) 

in agreement with ref. [9] up to an irrelevant overall phase^. 



(42) 



^ After correcting the sign in the second term of eqn. (3.18) of ref. [9]. 
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We have verified that eq. (37) has the correct colhnear hmits, and have also verified it 
numerically against a set of light-cone recurrence relations up to n = 12. 



V. THE A{ls,2+,...,m-,{m + l)+,...,ns) AMPLITUDE 

In this section, we will generalize the result from the previous section, and obtain an 
expression for the amplitude with a negative-helicity gluon not color- adjacent to one of the 
massive scalar legs. Here, we will use an (m, m — 1) shift, 

\m~) \m~) — zKm — 1)~) , 

|(m - 1)+) ^ |(m - 1)+) + z |m+) , (43) 

Our target amplitude is then given by the following recursive form 
^„(1„ 2+, . . . , m-, (m + 1)+, . . . , (n - 1)+, n,) 

= Amih, 2+, . . . , (r^l)+, -Klim-i),s)l7-2 

X ^„_^+2(-^it.(^_i),,, (m + 1)+, ...,ns) 



n-l 



+ ^n+m-i(ls) 2"*", . . . , (m — 1)"*", (i -|- 1)"*", . . . , rig)- 



=m+l (i+l)---(m-l) 
+As{{m - 2)+, i^iy, -k^^_^^^„^_^^)-^ 

^{m-2){m-l) 

X A„_i(l„ 2+, . . . , (m - 3) + , i^J„-2)(m-l)' • • • ' 

= C/i(l„ 2+, . . . , m-, . . . , (n - 1)+, n,) + ^3((m - 2)+, (^1)+, -K^^^2)im-i)) 

X -JH-^ ^n-i(ls, 2+, . . . , (m - 3)+, , m", . . . , n+) (44) 

-"■(m-2)(m-l) 

where as indicated Ui is defined as the first two terms. It depends only on gluonic ampli- 
tudes, all-plus massive-scalar amplitudes and m = 1 lone-negative-helicity massive-scalar 
amplitudes. We obtained all-multiplicity solutions for these in previous sections, and from 
them we can obtain the all-n form of this function. The Ui term corresponds to the h = 2 
terms in eq. (47). 

For the last term in eq. (44) we can, as before, iterate the recurrence relation eq. (44), 
whereupon we obtain 

Anils, 2+, . . . , m", (m + 1)+, . . . , (n - 1)+, n,) 
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^2^a_3(-i^[;],(m-r)+,ie+_,,) 

j=l r=2 



(m— r)---(m— 1) 



xC/i(l+ 2+,...,{m-j- 1)+, K+, m-,..., n,), (45) 

where K[r] and K^^] are as defined in eq. (32), and where we have used 

An-m+ai^s, K^^^^,rh~ , . . . ,ns) = C/i(l„ X[+ _2], m~, . . . , n,). (46) 

The amphtude we seek is now defined entirely in terms of Ui. Using this iterated recur- 
rence form of the amphtude along with the explicit form of Ui and also eq. (36) we arrive 
at our final form for the amplitude, 

An{ls, 2+, . . . , m", . . . , (n - 1)+, n^) = 



m—l n—2 



((2---(n-l))) 
/o(n,m;Zi)7o(n,m;Z2) 



S £ ^^"^"1 ^1— + ^2...mh \{l2 + 1) + ) 



+ 



(m— ii + l)---n (m-/i+l)-(n-l) l("^ ~ 

1 

m---(;2+l) (■^('2+2)---(r^-l)-^(rrjr-ii+l)---n ~l~ ^n-^(Trjr-ii+l)---(r!r-l))-^^(rrjr-ii+l)---m 1^ 

/i(n, m, Zi, ^2 + 1; h)f2(n, m, h, k + 1) (m"! t m-{h+i)t {m-h+i)-m 1™+)^ 



X 



X 



(47) 



-f^(^„i-ii+l)...(i2+l) ((^2 + 1) I $^ {rn-h+l)-{l2+l)^ {m-h+l)-m 1"^ 
1 



+\ 



X 



(m I ^„...(,2+i)^(rr»-ii+i)...(;2+i) |(m - Zi + 1)^ 
which is valid for 2 < m < n — 1. (For m — n — 1, one should use either eq. (37) or (38); for 
m = 2, one should use those formulae after refiection.) In the above equation. 



fo{n,m;i) 



fo{n,m]i) = 
and 

fi{n,m,li,l2;i) = < 



((m -i){m-i + 1)) G{2,m- i + l;m)/ L^rn-i+iy-n 

1, 

((i + 1) (i + 2)) C;(n; i + 1, n - 1; m) 

((« + 2)" I t2-{i+\)tx...m + ht2-m 1^^) 

' G'o(2,m-/i + l;/2,n-l;2;m) 
((m — -i) (m — i + 1)) Gi(2, m — i + 1; I2, n — 1; 2; m) 



i = m — 1; 
-, 2 < i < m - 1. 

(48) 

i = n — 2; 
-, m < i < n — 2. 



Li2{{m-i) |^(„_i+i).../2^^...z2 



i = m — 1; 
, 2 < i < m - 1. 
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(49) 

t ^ 1 1 .A\ — ) ("^~|-^m-{n-l)^n-^(m-ii+l)-{n-l)-?^(m-/i+l)-ml"^'*") 

1 „ ^ -<'(' + ^'> — , ™ + i<i<„-i. 

The attentive reader will notice that the dimensionality of the expressions for /i and /2 is 
different in each of the two cases. Nonetheless, all terms in eq. (47) have the same dimension, 
because this difference is compensated by the inapplicability of the first replacement in 

eq. (54) when I2 = n — 2. 
We have also defined 

Gu{ai,bi;a2,b2;c;m) = 



X 



(Wl+I ^ai-(«,i-l)^1^2-m , U = 0, 

{wi^\$a^-{wi-i)h\ar) , u = l, 



Wj =62 



where the sum for each Wi is over a set of momenta 5",, defined as follows. Define the set 
'S'o = {oi, ■ ■ ■ ,bi,a2, ■ ■ ■ , 62}, let succ(a) be the element following a, and pred(a) the element 
preceeding a in S'o. The set Si is then defined by omitting the first and last two elements of 
-S'o, 

5*1 = {succ(ai), . . . , 61, a2, . . . , pred(pred(62))}- (51) 

Each subsequent set Si depends on the value of Wi^i, it will contain all elements in Sq — 
{oi, . . . , 61, 02, . . . , ^2} following the element after Wj-i, and through two elements prior to 
62- This is the generalization of the sum for Wi starting at Wi^i + 2 (and ending at Wj+i — 2). 
That is. 

Si = {succ(succ('u;i_i)), . . . , pred(pred(62))}- (52) 

Note, however, that the sums of momenta Ki...j are over consecutive momenta, not restricted 
to So- 

The prime signifies that in addition, we must make the following replacements whenever 
Wk is inside a bra or a ket. For Wk = bi, 

(^fc |-^(«)fc+l)-(«;fc+i-l) ("^"|-^m-a2^bi-a2^«;fc-(«;fc+i-l) 



m---a2 
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and for Wk = a2, 



t {wk-i+l)-{wk-l)\'^k) t{wu-i+\)-wufm-a2\T^'^) (54) 

These replacements should also be applied to {wi \ in the last factor. That is if wi = 

(Vl^ai-K-l) {m'lth^-mtar-iwr-l)-, (55) 

and if Wi = 02, 

(W^l^ai-(«;i-l) -(^"l^m-aa^ai-m- (56) 

This definition of G is a generalization of that given in eq. (39), the two are related 
via G(ai, 61; m) = Gi(ai, hi\bi,bi\\\m) as {ai, . . . , 61} |J{^i} = {'^i; ■ ■ 1 ^1} only the 
second replacement from eq. (53) is relevant. Similarly G is defined as 

L(62-a2+l)/2j ba-l ^^'^ I -\lt I - \ 

G(«;a„6.;m)= ^ (-m,y (J] ^^^^^^„^,^^^„...„ ) (5^) 

{■<"i}i=l=a2+2 ^ 
Wi=u'i+l+2 



X 



(tui+l)---(n— 1) L \h2-) 



{^'j\ = (m I -Km- ■•02 



There is only a single massless contribution in eq. (47); it arises from the {li — m — ljh — 
n — 2) term and is given by 

((2 ■ ■ ■ (n - 1))) {m-\ Jt„...(„_i)^„ |(n - 1)+) 

{m~\ ^ntm-(n-l) \m''f .^„. 

X / -, TT / -I ri^ (58) 



which is readily seen to reduce to the expected MHV amplitude in the massless limit. Using 
eq. (47) we can derive the form of the amplitude in the five point, m = 3 case 

^5(ls,2+,3-,4+,5,) 

iml [2 4]^ 



i(3- 




2-)'(3- 




4-)^ 


Li2 (2 3) (34) L45 (4+11^5^34 


2-) 



34 1^ / 

this result matches exactly that given in ref . [9] . For the six point case with m = 4 we have 

^(l„2+,3+,4-,5+,6,) = 

16 



i{A-\ 




23 


|4+)'(4-| 




(2 3) (3 4) (4 5)L56 (5+| 


4^2...5it23 |4+) (4-| 




|2+) 



im^ [2 3] (4- 


-^56 


|3")'(4- 


^"6 


15-)' 


Li2L4...6l>56 (4 5) (2-| 
[3 5]^( 


5+1 


^23 
^2.. 


|4+)(5+| 

•4^1 |2 ) 


^6-^45 


|3-) 



Li2 [34] [45] (2-| |5-) (5+| ^^45 13") A^^..^ 



im^ (4 




-^23 


15-)^ 


i^|...5i^|...4 (2 3) (3 4) (5+1 


l^6^2...5^23 |4+) (2-| 


^34 


|5-) 



(60) 



Although this result has a different form from that given in ref . [9] , we have checked that the 
two agree numerically up to an overall phase^. We have also compared the result eq. (47) 
numerically to amplitudes computed via light-cone recurrence relations, and found complete 
agreement through n — 12. 

VI. CONCLUSIONS 

On-shell recursion relations have emerged as a powerful method for deriving analytic 
expressions for tree-level amplitudes. In this paper, we have given all-multiplicity solutions 
to the BGKS recursion relation for amplitudes with a massive scalar pair for the two simplest 
helicity configurations. These are the amplitudes with an arbitrary number of positive- 
helicity gluons, and either no or one gluon of negative helicity. (Of course, the corresponding 
amplitudes with an arbitrary number of negative-helicity gluons, and up to one gluon of 
positive helicity, can be obtained by spinor conjugation.) We have checked the principal 
results, eqs. (12), (37), and (47), against a numerical implementation of the older light-cone 
recursion relations, through n = 12. 

The primary application of these results is to computations of one-loop QCD amplitudes 
within the unitarity-based method. The massive scalar amplitudes are equivalent to ampli- 
tudes with the scalar legs computed fully in D dimensions. These in turn can be used to 
obtain rational terms in amplitudes, which have no cuts in four dimensions, but do have 
cuts in D dimensions. They should also be useful for proving general factorization proper- 
ties in complex momenta upon which the recursion-relations part of the unitarity-bootstrap 
combined approach relies [30]. For such applications, a relatively compact analytic form is 
desirable, and the recursion relations allow one to obtain such a form. 

^ After correcting a typographical error in eqn. (3.21) of ref. [9]. 
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While the all-plus amplitude did require guesswork for an all-n ansatz, it is interesting 
to note that the amplitudes with a ncgativc-hclicity gluon could then be derived without 
having to guess an ansatz for the general form. This would not be the case using the older, 
off-shell recursion relations [27, 28] . While the over-all complexity of the amphtudes would 
increase with increasing number of negative-helicity gluons, the same approach of unwinding 
three-point vertices could be used in such cases as well. These too would not require new 
ansatze for massive-scalar amplitudes. The calculations in this paper also illustrate the 
extent to which calculations beyond those in purely massless gauge theories are feasible. We 
would expect, for example, that amplitudes with massive quarks or vector bosons instead of 
massive scalars would yield to a calculation of similar nature. 
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APPENDIX A: COMPUTATIONAL COMPLEXITY OF BERENDS-GIELE RE- 
CURSION RELATIONS 

In order to use QCD amplitudes in evaluating experimental data, one must evaluate them 
numerically. How many operations are required to do so? Even at tree level, each ampli- 
tude, and indeed each different helicity amplitude, require a factorial number of Feynman 
diagrams. This might seem to imply that a factorial number of operations are required; as 
we shall show here, a better organization of the calculation can dramatically reduce this 
complexity. 

The complete squared amplitude can be written in terms of color-ordered amplitudes. 
At leading order in the number of colors, which is all we will consider here, computing 
the amplitude does require a sum over a factorial number of color orderings. However, 
for integrated quantities, one can use the symmetry of final-particle phase space to reduce 
this to a sum over a linear number of orderings. We are thus left to consider the cost of 
computing a given color-ordered amplitude. 

There are 0(2"/^) independent helicity amphtudes, so clearly the complexity of computing 
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all of them must scale at least exponentially in the number of external legs. There are an 
exponential number of color-ordered Feynman diagrams for any given hclicity amplitude, 
which would suggest that we still require an exponential number of operations. 

In this appendix, we show that the number of operations required to evaluate a typical 
helicity amplitude is only polynomial in the number of external legs, so long as the calculation 
is organized properly. To do so, we will make use of a light-cone version [28, 33] of the 
Berends-Giele recursion relations [27]. (The complexity is polynomial for any version of 
this type of off-shell recursion relation, suitably organized; only the prefactors differ.) In 
our estimates, we take each basic operation — addition, multiplication, division — to have 
constant complexity, that is 0{1). 

Let J^(l, . . . ,n) denote the (amputated) current with n on-shell external legs, and one 
off-shell external leg of the given helicity (defined using a hght-cone gauge vector q). We 
can write a recursion relation for this current, 

J- (l,...,n) = 

-EE Kr^2 J--(l,...,Ji)J--(ji + l,...,n) 

ji=l CTi,2=± ^•••■'1 (ji+l)---n 

ji=i j2=ii+i o-i,2,3=± ^--^'i 0i+i)-i2 (i2+l)-n 

X J-'^^(ji + l,...,j2)J-"^(j2 + l,...,n) 

In order to compute an n-point amplitude. An, we compute J'^(l,...,n — 1) with the n-th 
leg on-shell instead of off-shell. This in turn requires the computation of currents with fewer 
on-shell external legs. We can organize this computation as follows: first compute all the 
required three-point currents; then the four-point currents, and so on through the desired 
n-point current. At each stage, only previously- computed lower-point currents appear, so 
we need compute only the vertices and propagators, then multiply factors and perform the 
sums. In practice, it is probably more convenient to use caching rather than precomputation. 
That is, use the recursive formula, but during the recursive descent, record the (numerical) 
value of each newly-computed current. The next time that current is required, use the 
previously-computed value instead of computing it anew. 

For each j-point current contributing to the final n-point current, the computation is 
clearly dominated by the four-point vertices and associated double sums. We can also 
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precompute all momentum sums that appear in our computation. These are always sums of 
consecutive momenta, so there are (!?(n^) different ones, and this precomputation is of 0{n^) 
in operations. We can also precompute required Lorentz and spinor products, which requires 
0{n'^) operations. Each four-point vertex then takes a constant number of operations to 
compute. There are 0{j^) different terms, each requiring a constant number of operations, 
and so the computation of each j-point current requires 0{j^) additional operations. 

For a given j, we need n — j + 1 different j-point currents: J(l, . . . ,j — 1) through 
J{n — j + 1, n — 1). The overall computational complexity is thus of order, 

- j + l)f ^ 0{n% (A2) 

which is indeed polynomial. 

This is possible because we can re-use information: the computation of Jn, for example, 
requires knowledge of both Jio(l, • • • , 10) and J5(l, . . . , 5). The computation of Jio(l, • • • , 10) 
in turn requires knowledge of J5(l, . . . , 5), but this five-point current is the same as the one 
appearing in the top-level sum for Jn, and thus need be computed only once. This saving 
would not be possible if we insisted on writing out an analytic expression for in terms 
of spinor products, because the expression for J5 would then appear (and effectively be 
computed) many times. 

We can obtain a more precise count of the number of operations as follows. Denote the 
number of operations required to compute a three-vertex V3 by C3, and that for a four- 
vertex V4 by C4. If all helicity arguments to a three- vertex are identical, it vanishes, so for 
a given choice of the current's helicity a, there are only three different helicity choices we 
must compute. A straightforward organization of the sums requires four multiplications or 
divisions and one addition for each one, for example, 

KfKi ^ KlKl ^ KfKi ' ^ ' 

where Ji = J(l, . . . , ji), Ki = Ki j^, etc. It is possible to reduce this slightly by combining 
terms as follows, 

-^[\/-+- Jf J+ + J^(V-^Jt + V-++Jr)], (A4) 

giving instead a count of 3c3 + 10 for a typical term, or (3c3 + 10) (j — 2) for a j-point current. 
(This is still an overestimate of the optimal count, because wc ignore the reduction due to 
terms containing a two-point current as it depends on the choice of external helicities.) 
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Only four- vertices with two negative-helicity and two positive-helicity arguments are non- 
vanishing, so here we again have three different hehcity configurations for each choice of the 
current's hehcity. This would give six multiplications and one addition for each term in the 
double sum over ji and j2 and 8(04 -|- 7){j — 2){j — 3)/2 overall. Here, we can reorganize not 
only each term but the double sum as well. 



j-3 j-2 

^ ^ [ KfKiKi ' KfKlKl ' KfKiKi 

jl=lj2=ji+l 1 2 3 1 2 3 1 2 3 



n 



i^2j^2 

J2=jl + 1 2 ^ j2=jl + l 2 ^ 



(A5) 



giving an overall operation count of 



^(c4 + 4)(j-2)(j-3) + 4(j-3) (A6) 

(again ignoring the reduction due to two-point currents). 
Overall, we then obtain an operation count of 

J2in - J + 1) [(3C3 + 10)(j - 2) + -(C4 + 4)(j - 2)(j - 3) + 4(j - 3) 

.=3 ^ ^ (A7) 

n(n- l)(n-2)(n-3) n(n-l)(n-2) (n - l)(n - 2)(n + 3)(3n - 4) 

= 8 + 2 + 6 

In contrast to the off-shell recursion relations analyzed here, in on-shell recursion relations, 
the reference momenta will change as the recursion descends. It thus seems likely that 
the computational complexity will still be exponential for a typical color-ordered helicity 
amplitude. Of course, for special helicity configurations or for a moderate number of external 
legs, the new representations may provide a more efficient computational method. 
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